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The light scalar meson spectrum is studied using the improved ladder QCD with the Ua{1) break- 
ing Kobayashi-Maskawa-'t Hooft interaction by solving the Schwinger-Dyson and Bethe-Salpeter 
equations. The dynamically generated momentum-dependent quark mass is large enough in the low 
momentum region to give rise to the spontaneous breaking of chiral symmetry. Due to the large 
dynamical quark mass, the scalar mesons become the qq bound states. Since the parameters have 
^~») ' been all fixed to reproduce the light pseudoscalar meson masses and the decay constant, there is no 

I free parameter in the calculation of the scalar mesons. We obtain = 667 MeV, Mag = 942 MeV 

, and Mfg = 1336 MeV. They are in good agreement with the observed masses of (t(600), ao(980) 

' and /o(1370), respectively. We therefore conclude that these states are the members of the light 

' scalar meson nonet. The mass of Kq is obtained between that of ao and fo and the corresponding 

state is not observed experimentally. We also find that the strangeness content in the a meson is 
about 5%. 
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I. INTRODUCTION 
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I About three decades ago, quantum chromodynamics (QCD) was conceived as the microscopic theory of the strong 
■ interactions. Since then, many aspects of QCD have been studied and QCD is established as the fundamental theory 
of the strong interactions. Because of the nonperturbative nature of the low-energy QCD, understanding the low- 
lying hadron structures from the viewpoint of the quark and gluon degree of freedom is one of the most challenging 
' problems. 

The lightest excitation on the QCD vacuum is the pion which is considered as a quark and an antiquark bound 
state in the pseudoscalar channel. Its mass (about 140 MeV) is off-scale light compared with other hadrons such as 
Qh! the /9- meson {nip ^ 770 MeV), the spin- flip partner of the pion, and the nucleon {rriN ^ 940 MeV), the three-quark 
^ I bound state. This can be understood by recognizing that the chiral symmetry is spontaneously broken in the QCD 
vacuum and the pion is the Nambu-Goldstone (NG) boson associated with the dynamical chiral symmetry breaking 
!>■ ' (DCSB). The DCSB is, therefore, among the most important aspects of low-energy hadron physics. It is believed to 
' be responsible for a large part of the constituent quark masses, which are introduced in the many constituent quark 
. models. 

If QCD lagrangian has no explicit chiral-symmetry breaking term, the mass of the NG boson associated with the 
DCSB should be zero. In order to explain the observed mass of the pion, one needs small explicit chiral-symmetry 
breaking terms, namely, current quark mass terms. The NG boson nature of the pion was first studied by Nambu and 
Jona-Lasinio (NJL) using the schematic model with the four-fermion interaction .1]. Since then, low-energy hadron 
properties have been widely studied using the NJL-type models On the other hand, the effects of the explicit 
breaking of the chiral symmetry on the pion properties have been systematically studied using the effective lagrangian 
composed of the pion field. This approach is called the chiral perturbation theory (ChPT) The success of ChPT 
approach supports the importance of the DCSB in low-energy QCD. 

When we start looking at the strange quark sector, we meet another problem. If one assumes that the up, down 
and strange quark masses are small compared with the scale of the DCSB, the number of the NG bosons is equal to 
the dimension of the coset space Ul{3) x /7i{(3)//7v'(3), namely, nine. The ninth heavier pseudoscalar meson is rj' and 
its mass is much heavier than the other octet pseudoscalar mesons. Weinberg showed that the mass of rj' should be 
less than \/3mT, if ?7a(1) symmetry were not explicitly broken 4j. Thus the Ua{^) symmetry must be broken. The 
key step to solve this Ua{^) problem was to realize that there is an anomaly in the C/a(1) channel. Namely the Ua{^) 
symmetry in the classical theory, i.e., in the action, is broken by quantum effects. In the following year, 't Hooft 
pointed out the relation between Ua{^) anomaly and topological gluon configurations of QCD and showed that the 
interaction of light quarks and instantons breaks the C/a(1) symmetry^. He also showed that such an interaction can 
be represented by a local 2Nf quark vertex, which is antisymmetric under flavor exchanges, in the dilute instanton 
gas approximation. 
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The effects of the anomaly on the low-energy QCD have been extensively studied in the l/iV"c expansion 

approach |^. In the Nq — *■ oo limit, the 1/4(1) anomaly is turned off and then the 77 and 77' mesons become the ideal 
mixing states. The flavor component of the ry is ^(mm + dd) with mri{Nc ~^ 00) — m-^ and the ?/ meson becomes a 
pure as state with m^,{Nc 00) ~ 2m|- — — (687MeV)^ |3|. The higher oder effects of the l/Nc expansion give 
rise to the flavor mixing between the ry and rj' mesons and push up the ry' mass. They were further discussed in the 
context of the ChPT [3] and the reasonable description of the nonet pseudoscalar mesons was obtained. 

The breaking 27V/ quark determinant interaction was introduced to the low-energy effective quark models 

of QCD. The low-lying meson properties have been studied 0, 0, O, 0, 0| using the three-flavor version of the 
NJL model with the Kobayashi-Maskawa-'t Hooft (KMT) determinant interaction The radiative decays of 

the rj meson have been studied in this approach Il5| and found that these decay widths are reproduced when the 
C/A(l)-breaking interaction is much stronger than the previous studies 9, idUl [12. IT^j l . It is further argued that the 
f/A(l)-breaking interaction gives rise to the spin-spin forces, which are important for light baryons |lft [l7l [l8l| . 

The dynamics of instantons in the multi-instanton vacuum has been studied by many authors, either in the models 
or in the lattice QCD approach, and the widely accepted picture is that the QCD vacuum consists of small instantons 
of the size about 1/3 fm with the density of 1 instanton (or anti-instanton) per fm^ p^ . In the instanton liquid model, 
the instanton plays a crucial role in understanding not only the J/a(1) anomaly but also the spontaneous breaking of 
chiral symmetry itself. 

Recently, the low-lying scalar mesons, J"^ = O"*", attracts a lot of attention by two reaso ns. l2d | (1) Experimental 
evidence for a {I = 0) scalar meson of mass around 600 MeV is overwhelming pH 123. l2.ll l24|. Especially the 
decays of heavy mesons show clear peaks in the tttt invariant mass spectrum. Including this a rather light isoscalar 
state, the light scalar mesons show strange mass patterns, i.e., cr(600) — ao(980) — Kq{?) — /o(980), where Kq is not 
confirmed. This pattern cannot be explained as a naive qq nonet, because the / = 1 ao states are almost degenerate 
with the second I ~ state /o, while the first / = state a is far below them. (2) The roles of the scalar mesons 
in chiral symmetry have been stressed in the context of high temperature and/or density hadronic matter [2^. It 
is believed that chiral symmetry is restored in the QCD ground state at high temperature (and/or baryon density). 
Above the critical temperature, ~ 170 MeV, the world is nearly chiral symmetric and we expect that hadrons belong 
to irreducible representations of chiral symmetry, if we neglect small mixing due to finite quark mass. The pion is 
not any more a Nambu-Goldstone boson, and has a finite mass and should be degenerate with a scalar meson, i.e., 
sigma. Another scenario was proposed rec ently . It is the vector manifestation where the rho meson becomes massless 
degenerate with pion as the chiral partner [2^. In order to make the situation clear, we consider that the light scalar 
mesons are the key particles. 

We have studied the effects of the C/a(1) breaking interaction on the low-lying nonet scalar mesons using the 
extended NJL model, in which the three-flavor NJL model is supplemented with the KMT determinant interaction 
[27^ . Why is the J7a(1) expected to be important in the scalar mesons? It is because the KMT interaction selects out 
the scalar sector as well as the pseudoscalar mesons and therefore the OZI rule may be broken significantly also in the 
scalar mesons ,28J. We have found that the C/a(1) breaking interaction gives rise to about 150 MeV mass difference 
between the a and qq mesons. We have also found that the strangeness content in the a meson is about 15%. The 
calculated mass of the 7 = 1/2 state (Kg) was about 200 MeV heavier than that of the 1=1 state (oq). 

The physics of the light scalar mesons seems to be directly related to the mechanism of the dynamical chiral 
symmetry breaking (DCSB). In the NJL model, the strong four-quark interaction gives rise to the DCSB and it leads 
the simple mass relation between the scalar meson mass (ms) and the dynamical quark mass (Mq) in the mean-field 
approximation, i.e., ms = 2Mq in the chiral limit. In the case of the instanton liquid model, the DCSB is caused 
by the instanton induced interaction. The scalar meson masses are rather sensitive to the instantion-antiinstanton 
interaction. The results in the fully interacting instanton ensemble are TOo- — 0.58 GeV and uiag — 2.05 GeV JJ]. 

In contrast with the instanton liquid model, the study of the QCD Schwinger-Dyson (SD) equation for the quark 
propagator in the improved ladder approximation (ILA) has shown that the spontaneous breaking of the chiral 
symmetry is explained by simply extrapolating the running coupling constant from the perturbative high-energy 
region to the low-energy region [23. Then, the Bethe-Salpeter (BS) equation for the J^^ = qq channel has 
been solved in the ILA and the existence of the Nambu-Goldstone pion has been confirmed [s^,!^. The numerical 
predictions of the pion decay constant and the quark condensate {tl>tp) are rather good. It has been also shown 
that the BS amplitude shows the correct asymptotic behavior as predicted by the operator product expansion (OPE) 
in QCD ,32] . The masses and decay constants for the lowest lying scalar, vector and axial- vector mesons have been 
evaluated by calculating the two point correlation functions for the composite operators ipMijj. The obtained values 
are in reasonable agreement with the observed ones |33| . 

Then, the current quark mass term has been introduced in the studies of the BS amplitudes in the ILA |34| and the 
reasonable values of the pion mass, the pion decay constant and the quark condensate have been obtained. It has been 
also shown that the pion mass square and the pion decay constant are almost proportional to the current quark mass 
up to the strange quark mass region. The effect of the Ua{^) anomaly is further introduced in the ILA approach by 
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the KMT interaction. The instanton size effects are taken into account by the form factor of the interaction vertices. 
It guarantees the right asymptotic behavior of the solutions of the SD and BS equations. The properties of the rj 
and rj' mesons have been studied by solving the coupled channel BS equations and the reasonable values of mjr, m^, 
n^ri' , /tt and {qq)B. with a relatively weak flavor mixing interaction (KMT interaction), for which the chiral symmetry 
breaking is dominantly induced by the soft-gluon exchange interaction [35l | . 

The purpose of this paper is to study the properties of the light scalar meson nonet in the improved ladder ap- 
proximation (ILA) of QCD with the breaking Kobayashi-Maskawa-'t Hooft (KMT) 6-quark flavor determinant 
interaction. In this approach, the mechanism of the dynamical chiral symmetry breaking is different from the NJL 
model and instanton liquid model. It has been shown that the Wilsonian non-perturbative renormalization group 
equation gives the identical effective fermion mass with that obtained by solving the Schwinger-Dyson equation in the 
improved ladder approximation [s^. We hope that the present study may shed light on the mechanism of the DCSB 
in the low-energy QCD. It should be noted here that the parameters have been all fixed in the pseudoscalar meson 
sector and therefore there is no free parameter in the present study. 

There have been many studies of the pion BS amplitude using the effective models of QCD and /or the approximation 
schemes of QCD [s^ lSSj . The main differences among these studies are the form of the gluon propagator used in the 
SD and BS equations. Of course, the behavior of the gluon propagator in the asymptotic region is well established and 
there is no differences. However, that in the infrared region is not well known. The simple infrared cutoff is introduced 
in the ILA, while in the approach given in 's^, the specific form of the infrared gluon propagator is assumed. 

Recently, attempts have been made to evaluate enhancement of the gluon propagator (as well as the vertex factor) 
from quenched lattice QCD data of the quark propagator 39, 40]. They found that the gluon propagator is required 
to have strong enhancement in the soft momentum region so that chiral symmetry is dynamically broken. Although 
the enhancement patterns vary among the models and parametrizations, they tend to show similar behavior when 
the solution of the SD equation, i.e., the quark mass function, is concerned. Thus we expect that the results for the 
mesonic bound state in the BS equation with the effective quark mass function qualitatively agree among the models. 

The paper is organized as follows. In Sec. ^ we explain the Lagrangian we use in the present study and derive the 
SD equation for the quark propagator and the BS equation for the scalar meson. Sec. Illll is devoted to the numerical 
results. Finally, summary and concluding remarks are given in Sec. IIVI 



In this section, we present the formulation of the improved ladder QCD with KMT interaction. The rainbow 
approximation is applied to the SD equation for the quark propagator and to the BS equation for the pseudoscalar 
and scalar mesons. Since the derivations of the SD equation and the BS equation for the pseudoscalar mesons have 
been given in the ref. [35l |. here, we present only the results. 



The improved ladder QCD is based on the ladder approximation which is improved by replacing the coupling 
constant by the running coupling constant. We employ the rainbow approximation of the SD equation for the quark 
propagator and the ladder approximation of the BS equation for the quark-antiquark bound states. Improvement is 
made by the use of the running coupling constant according to the Higashijima-Miransky [2^l30l l33| method. 

Under this approximation, the gluon exchange interaction £ge becomes 



II. FORMULATION 



A. Improved ladder QCD with KMT interaction 
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FIG. 1: The Schwinger-Dyson equation 
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In eq.|0J the infrared cut-off tip is introduced. Above tip, g'^iPs) devefops according to the one loop solution of the 
QCD renormalization group equation, while below to, g^(p%) is kept constant. These two regions are connected by 
a quadratic polynomial so that g'^{p%) becomes a smooth function. Here Nc and Nf are the number of colors and 
active flavors respectively. We use Nc — Nf = 3 in this study. 
The KMT interaction £kmt is given by 
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KMT 



Qj^gfif2 fa ^919293 



<P^+P-+P^+^^+i-+i^^-wipi,p2,P3,qi,q2,q3) 



X { [^31 (gi)] [^92 (P2)V'/2 (92)] [^93 (P3)V'/3(a3)] 

+3 [^91 (91)] ^2 {P2)ibi^h (92)] [^93 (pshbi^h (93)] } 

where fi,gi, ■ ■ ■ are flavor indices, e denotes the antisymmetric tensor with e""^" 
symmetry and also mixings quark flavors. We introduce a weight function 



w{pi 



,93) = exp(-K(pi + h gl)) 



(9) 

1 This interaction breaks the J7a(1) 



(10) 



so that the KMT interaction is turned off at large momentum region. Then the asymptotic behavior of the ILA are 
kept consistent with the perturbative QCD. The parameter k is taken as k = 0.7[GeV~^]. This value corresponds to 
the form factor of the instanton of the average size p ~ 1/3 [fm] |35l |. 



B. Schwinger-Dyson Equation 

The SD equation in the momentum space is written as 

iSF-\q) - tSo-\q) = ~Cf I f{-q\-p^)lD^-{p-q)-f,SF{p)-1, 

Jp 

-Gn I wi-q'-p' - A:2)tr(DC)[5^(p)]tr(DC)[5^(fc)]. 



(11) 
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FIG. 2: q% dependencies of the mass function of the light (q) quark, Bg, and the stranger (s) quark, Bg, for Go ~ and 
75[GeV-^] 



where tr^^^^ means the trace of the Dirac and color components and Cp = {N^ — l)/2Nc- Here, the flavor antisym- 
metrization is assumed in the second term of the RHS. This equation is shown diagrammatically in Fig. ^ 
Generally the quark propagator is parametrized as 



(12) 



In the Landau gauge, it can be shown that the solution satisfies A{—q'^) — 1. Then the SD equation becomes the 
integral equation only of the mass function B{q^). Our numerical results of the mass function are shown in Fig(21 
The quark masses are renormalized as 



niq — Z^^^niqj^, 



where we take the renormalization condition as 



= 1 



nisfi—O 



We define the quark condensate as follows 

imm) = - [ tiis^ip)] 
jp 

where 
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FIG. 3: The Bethe-Salpeter (BS) equation 



C. Bethe-Salpeter eqation 



To treat the scalar and pseudoscalar mesons as quark-antiquark bound states, we use the homogeneous Bethe- 
Salpeter equation. The BS equation in momentum space becomes 

Sp\q+)x{q;P)Sp\q-) - -zCp [ f{-q\-e)lD^^-{q~kh^x''{k;Ph. 

Jk 

~2^Gn jjw V - - - ^) trf^^) [S^{p)] 



{75tr(°c) p^)] ^ t^(DC) 



(19) 



where q± — q± Pb) denotes the BS amplitude and Pb is the momentum of the meson. Like the SD equation, 

the Gd term is antisymmetrized in the flavor space. All the momentum integrals are regularized by the cut-off Ayvi 
which is chosen as Auv = 100[GeV] so that solutions of the BS equation do not depend on it. This equation is shown 
diagrammatically in Fig. |21 

For the scalar meson, the BS amplitude can be parametrized in terms of four Dirac structures with appropriate 
flavor structures, 
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X^(fc;P) - lcY['l'sik;P) + rp{k;P)^ + rQ{k;P)If + <p^^{k-P)-m - ^ 
The normalization condition of the BS amplitude is given by 



(20) 



d^q d 

* / 7Tr^x«i™i('?;^')xm2n2(9;^')-P''7^ {'S'F;:2\,i(g+)5'F;;^,„.,((7-)} = -2P^ (21) 



dPt^ 



where the indices toi, ni, . . . are combined indices in the color, flavor and Dirac spaces. 

In the numerical computation, we solve the BS equation in the Euclidean momentum region. Then the physical 
solution, which corresponds to negative P^, is obtained by extrapolation from the P| > region. It can be done in 
the following way. First, we rewrite the Euclidean BS equation in the form 



bAiq]PE)-- / MABiqE]kE;PE)(l)Bik;PE) 

JkE 



(22) 



where (I)a and (pB denotes a set of amplitude. This equation should not have a solution at Pj. > because, if there is 
one, it is a tachyon solution. Therefore we instead solve an eigenvalue equation 



X{P^)<j)A{qE;PE)= [ MAB{qE;kE]PE)MkE]PE) 

JkE 



(23) 



for a fixed P| > 0. Then we extrapolate the eigenvalue A to P| < as a function of Pj. and search for the on-shcU 
point A(— M^) = 1. We employ the quadratic function of P| for the extrapolation. 



III. NUMERICAL RESULTS 



A. Parameters 



In the present approach, there are seven input parameters: the bare quark mass ruq^ for the up and down quarks, in 
which we assume isospin symmetry m„ = m^, the current quark mass TOsq for the strange quark, the scale parameter 
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TABLE I: The numerical results of the BS equation. 



mass (MeV) 


TT 


136 


A = 95 [MeV] 


V 


515 


mixing angle —20.0° 


v' 

K 


982 
517 


mixing angle —26.2° 


a 


667 


mixing angle —68° 


ao 


942 




/o 


1336 


mixing angle —83.9° 


7^0* 


1115 





TABLE II: The values of the parameters and the obtained observable quantities. 



m,^(2[GeV])[MeV] 


m,fl(2[GeV])[MeV] 


AQCD[MeV] 




U 




fi:[GeV"^] 


4.5 


150 


600 


0.204 


-6.89 


75 


0.7 



of QCD running coupling constant Aqcd, the infrared cutoff iip, the smoothness parameter to, the strength parameter 
of the KMT Gd and the parameter of the weight function of the KMT k. We choose these parameters to reproduce 
the observables of the pseudoscalar mesons and then we apply them to the scalar mesons. 

The quark masses JTigoi 'tIso are chosen so that the renormalized masses at the momentum scale /i = 2[GeV] become 
the m-q^ — 4.5[MeV] and ra^ji = 150[MeV], respectively. The k parameter is taken as k = 0.7[GeV~^], which 
corresponds to the instanton of the average size 1/3 [fm]. The other parameters Aqcd, ^if, io, Gd are chosen as the 
pseudoscalar meson masses M^-, M^, M^i and the pion decay constant are fitted to the experimental values. The 
parameters that we use are Aqcd — 600[MeV], to = —6.89, tip — 0.204 and Gd — 75[GeV~^]. These parameters give 
= 136[MeV], = 515[MeV], M,,/ = 982[MeV] and = 95[MeV]. These are in agreement with experimental 
values in less than 6% of deviation. Table HTl summarizes all the values of the parameters. 

Although Aqcd is somewhat larger than the standard value Aqcd = 100 ~ 300 [MeV], a large Aqcd is necessary in 
the ILA to generate the desired dynamical chiral symmetry breaking (D^SB). It is known that by using the running 
coupling constant obtained by the higher loop calculation, the smaller Aqcd gives rise to the correct size of the DCSB. 
On the other hand, the Wilsonian non-perturbative renormalization group analysis has shown that the effects beyond 
the ladder approximation also reduce Aqcd j4l| . 



B. Solution of the SD equation 

The numerical solutions of the SD equation are shown in Fig|21 The values of the quark condensates and mass 
function at P| = are given in the Table UTII 



C. Solution of the BS equation for the Scalar mesons 

As mentioned above, since the parameters of this approach have been chosen using the observables of the pseu- 
doscalar mesons, the following numerical results of the scalar mesons are parameter free predictions. 

Let us now start the discussion of the solutions of the BS equation. Our numerical results for the scalar mesons are 
summarized in Table ^ The dependence of the mass spectra on the strength of the KMT interaction is shown in the 
FigH 

First, the dependencies of the masses of a and on Gd look qualitatively same as the NJL results shown in [2^ . 
In the NJL calculation, the parameters are chosen so as to reproduce the M^r and at each Gd- In contrast, in the 



TABLE III: The values of the mass function at P\ = and the quark condensates. 





Bs(0) 




0.616[GeV] 


0.778[GeV] 


0.259[GeV] 
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FIG. 4: Dependence of the mass spectrum of the scalar and pseudoscalar meson nonet on the strength of the KMT interaction. 



ILA approach we change Gd independently. However, since and depend weakly on G^j, the results of ILA 
show similar behavior as those of NJL. We note that the mass of ag grows as Gu increases, while the tr mass is almost 
constant. 

The a meson mass is predicted as about 670 [MeV]. This rather small a meson mass is interesting. In the case of the 
NJL model, the a meson mass is determined to be close to twice of the dynamical quark mass. On the other hand, in 
the ILA approach the value of the mass function at q^^ = 0, Bg{q^ = 0) is about 616[MeV], which is comparable to the 
a meson mass. Recently, Maris calculated the a meson mass in the SD and BS approach with the infrared enhanced 
gluon propagator and the obtained mass is 671 [MeV]. Although there are such differences, the properties of the 
physical observables agree in these calculations. 

Concerning the oq meson, we obtain Ma^ = 942 [MeV]. This result is comparable to the experimental value 
984. 8± 1.4[MeV]. We obtain significant mass splitting between the a and qq, about 275[MeV]. We conclude that the 
observed a — qq mass splitting can be explained as the J7a(1) symmetry breaking effects. Recent study of the light 
flavor scalar mesons in the instanton liquid model showed similar result [4^ . 

The obtained mass of /q is 1336 MeV. We therefore identify this state with /o(1370) whose T- matrix pole position 
is (1200 - 1500) - «(150 - 250) MeV. The calculated mass of /o is about 400 MeV larger than that of gq. The mixing 
angle of /o is —83.9° and it means that this state is close to the flavor octet state. In such a case, the effects of the 
KMT term on /o state is almost same as that on qq and therefore the origin of the mass difference between /o and gq 
is considered to be mostly the symmetry breaking effects by the strange quark mass. 

The /o(980) state is observed between o-(600) and /o(1370). From the present study, it seems unlikely that /o(980) 
is a simple gg-bound state. Literatures have suggested that /o(980) may consist of g^g^ [ij, [i^, Eg] or may be a 
KK molecular state Those studies often indicate that the ao(980) is also a four-quark state. It conflicts with 
our picture. Further study is necessary to make the situation clearer. 

The dependences of the scalar meson spectrum on the strength of C/a(1) breaking interaction are understood as 
follows. When = 0, cr meson is ideal mixing state -^{uu + dd). As the Go becomes large, the a approaches 

the flavor singlet state {uu + dd + ss) . The a meson mass seems to become large because of the increase of the 
strange component. On the other hand, the KMT interaction is attractive for the flavor singlet state, consequently 
the a meson mass hardly change. In the case of ao, since the KMT interaction dose not induce the flavor mixing, oq 
is flxed on the "^(wu + dd) state. As Gd becomes large, the uq mass increase by the repulsive KMT interaction. As 

Gd becomes large, the /o changes from ss to -^[uu + dd ~ 2ss). The decreasing of the strange component and the 
repulsive KMT interaction compete. At the small Gd, they are balanced and the fo mass increases at the large Gd 
where the flavor of the fo is sufficiently mixed. 

There is a shortcoming in the present approach. The solution of the BS equation is obtained at the Euclidian 
momentum region. To obtain the physical mass and mixing angle, we have to extrapolate the solution from the 
space- like to the time-like region. This is carried out by extrapolating the eigenvalue function A(P|) in Eq. H23|) to 
negative until it hits A(P|.) = 1. The graphs of the extrapolation about the masses of the ao, cr, fo are shown in 
the FigOWe therefore anticipate moderate ambiguity in the extracted masses especially when the mass is large. 
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FIG. 5: Extrapolation of the eigenvalue A with respect to for ct, ao and /o mesons. The solid line shows the extrapolation 
of the (7 and the doted line shows the extrapolation of the ao and the dashed line shows the extrapolation of the /o. 
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FIG. 6: Depnendence of the mixing angle of the a and /o meson on the strength of the KMT interaction. 



Next, we consider the mixing angles. We introduce the matrix elements and which are defined by 

/A" 
(i''a;(0|g— (7(a;)|scalar meson state) 



tr 



A-^ 



(24) 
(25) 



Since these S values extract the particular flavor component of 05 which is the main component of the BS amplitude 
at the origin, we employ and 5° to determine the ratio of the octet and the singlet components. Accordingly we 
define the mixing angles of the scalar mesons as 



tan 6„ 



tan0 



/o 



'-'0 

OCT 



fo ■ 



(26) 
(27) 



The results are summarized in Tableland Fig. El 

As Gd increases, the mixing angle moves towards —90°, where a becomes the purely flavor singlet state: -^{uu + 

dd-\- ss). The obtained angle at Gd = 75[GeV~^] is —68.0° and is slightly smaller than the result of the NJL model, 
—77.3°. This angle corresponds to about 5% mixing of the strangeness component in a. 
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D. Solution of the BS equation for the Strange meson 

In this section, we discuss the strange scalar meson Kq. Here we employ an approximation in order to avoid the 
technical difficulty coming from the ambiguity in defining the center of mass coordinate as this meson consists of the 
a strange quark and a nonstrange quark. Instead of treating the asymmetric BS equation, we solve the symmetric BS 
equation for the quarks of mass — 77.25[MeV]. We apply this approximation to the kaon and obtained the 

reasonable kaon mass, mx = 494 — 498 MeV. 

The results are summarized in Table 13 and Fig. 01 The obtained mass is 1115 MeV, which is about 173 MeV 
larger than that of oq and about 221 MeV smaller than that of /q. If all the oq, Kq and /o states are assumed to 
be the flavor octet qq states, the simple quark model calculation predicts the Gell-Mann-Okubo mass formula. Our 
results deviate from both the linear mass formula, 3(M/(, — Mk') = Mk* — Ma^, and the quadratic mass formula, 
mi - Ml,) = Ml, -Ml. 

Unfortunately, the corresponding light 7=1/2 scalar meson is not observed. The observed mass of the /irg(1430) 
is 1412 ± 6 MeV and larger than our result of Mf^. It therefore seems to be rather difficult to identify i4rQ(1430) with 
our state. We do not consider ifQ(1430) as a member of the light scalar nonet states. We hope the re-analysis of 
the experimental data using the chiral effective model including the effects of the t/A(l)breaking interaction will shed 
light on the problem of the missing state in the 1 = 1/2 scalar channel. 

IV. SUMMARY AND CONCLUSIONS 

We have studied the spectrum of the light scalar nonet mesons using the improved ladder approximation (ILA) of 
QCD with the Kobayashi-Maskawa-'t Hooft (KMT) Ua{^) breaking interaction. We choose parameters to reproduce 
the masses and decay constants of the pseudoscalar nonet mesons and apply those to the scalar nonet mesons. 

The ILA of QCD is an approximation that is consistent with chiral symmetry and consists of the rainbow approx- 
imation of the Schwinger-Dyson equation and the ladder approximation of the Bethe-Salpeter equation. Using this 
approach, we analyze the scalar meson spectrum quantitatively. 

We have obtained M„ = 667 MeV, Mag = 942 MeV and Mjg = 1336 MeV. Considering the ambiguity due to the 
extrapolation from the Euclid momentum, they are in good agreement with the observed masses of (t(600), ao(980) 
and /o(1370), respectively. We therefore consider that these states are the members of the light scalar meson nonet. 
This identification is different from the conventional one. The key ingredient is the instanton-induced breaking 
interaction. It gives rise to the symmetry breaking and flavor mixing effects on the scalar mesons as well as the 
pseudoscalar mesons. 

We have obtained the strangeness content in the a meson of about 5%. This ss mixing may be tested, for instance, 
by analyzing the tt tt decays of heavy mesons carefully. 

The obtained Kq mass is 1115 MeV. The corresponding state is not observed. The observed i4rQ(1430) is heavier 
than our result of Mfg and therefore we do not include 7^0(1430) in the light scalar nonet. 

It should be noted that the scalar isoscalar state can be mixed with the scalar glueball state. Such effect is not 
taken into account here. There is another important effect to be discussed. The light scalar nonet can be generally 
coupled to the intermediate states composed of two pseudoscalar mesons rather strongly. If the light q'^(f states exist, 
this coupling effects should be very important. The extension of the Bethe-Salpeter equation to the q'^q'^ system seems 
to be rather difficult. One way is to apply the bilocal bosonisation technique [i^. Another direction of the further 
study is the extension to the finite temperature and/or density jH^I- In the present study, the mass of the 7=1/2 
scalar meson is not reproduced well. In order to make the structure of the light scalar mesons clearer, the analysis of 
the experimental data using the framework including the Ua{^) breaking interaction is necessary. 
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